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^ ■ Abstract 
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ON ■ In this paper we re-examine the one- dimensional interaction of electromag- 

J>~>' netic and ion acoustic waves in a plasma. Our model is similar to one solved 

V ■ 

Q ■ by Rao et al. (Phys. Fluids, 26, 2^88 (1983)) under a number of analytical 

^ : 

■ approximations. Here we perform a numerical investigation to examine the 

• stability of the model. We find that for slightly over dense plasmas, the prop- 

^ • agation of stable solitary modes can occur in an adiabatic regime where the 

ion acoustic electric field potential is enslaved to the electromagnetic field of 
a laser. But if the laser intensity or plasma density increases or the laser fre- 
quency decreases, the adiabatic regime loses stability via a transition to chaos. 
New asymptotic states are attained when the adiabatic regime no longer exists. 
In these new states, the plasma becomes rarefied, and the laser field tends to 
behave like a vacuum field. 
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I. INTRODUCTION 



The recent developments of new laser technologies allow the creation of strong pulses that 
can propagate in a plasma, either to accelerate particles or induce laser assisted fusion 
0. Now, if a laser pulse interacts with a plasma the possibility exists of nonlinear wave 
coupling involving the pulse itself and nonlinear plasma modes [0,^. Since the variety of 
these nonlinear plasma modes is large, and since each mode exhibits a considerable richness 
regarding amplitude, polarization and frequency range [^,0, the nonlinear mode coupling is 
a feature to be appreciated with care. 

One important example of nonlinear mode coupling concerns the interaction of rela- 
tivistically strong electromagnetic waves with Langmuir waves These analysis are 

restricted to simple wave solutions with well defined propagation velocities which are super- 
luminous in an unmagnetized plasma. For these superluminous waves it is shown that the 
wave dynamics can be viewed as similar to the coupled dynamics of two nonlinear oscillators 
with their natural frequencies given by real numbers. One of the oscillators describes the 
transverse electromagnetic field and the other the longitudinal electrostatic field. The dy- 
namics is found to be predominantly integrable. Some nonlinear resonant island chains are 
present in the appropriate Poincare maps p[ , but their overlap is so small that the resulting 
trajectories are mostly regular. 

In another important range of subluminous wave velocities, around the ion acoustic 
range, the laser wave can couple to ion acoustic plasma modes. In this case the coupled 
waves propagate with velocities close to the ion acoustic velocity which is much smaller than 
the velocity of light. Here the ionic dynamics plays a crucial role and the resulting dynamics, 
in principle, bears no resemblance to the dynamics described in the previous paragraph. A 
good deal of analytical work has been done in Ref. to show that if one assumes once more a 
common and constant propagation velocity, the system becomes again equivalent to a pair of 
nonlinearly coupled oscillators. A novel feature studied here is that when the laser frequency 
is smaller than the average plasma frequency, wave localization can occur such that one has 
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the formation of intense solitary pulses. In terms of the system of equivalent oscillators, 
the possibility of solitary pulses occurs when by varying the appropriate parameters the 
linear frequency of that oscillator describing the laser wave becomes imaginary. Using the 
language of nonlinear dynamics, this change in the character of the frequency occurs because 
the central elliptic point present in the appropriate phase space turns into a hyperbolic point 

The work by Rao et al. ^ utilizes powerful approximation techniques, but is essentially 
analytical as commented before. Therefore one would like to have some information on 
the stability of the solitary pulses thus formed, and this is what we do here. The stability 
issue has been already addressed several times over the past years, but focus has been 



preferentially directed upon the linear perspective ||TT|. However, if one wishes to have some 
additional insight on the nonlinear development of these instabilities, the most appropriate 
tools of investigation appear to be the Poincare maps mentioned above. With these maps 
one records the phase space coordinates of one of the oscillators, as one of the coordinates 
of the other crosses its zero with a definite sign for the derivative. Since the system is two 
degrees of freedom and since there exists a conserved Hamiltonian for the system, the point 
recorded on the map gives all the relevant information for the dynamics. Regular motion is 
associated with smooth curves of the Poincare maps and chaotic (or nonintegrable) motion 
is associated with an erratic distribution of points representing the trajectory. We find 
here that the stability of the solitary pulses is quite limited. In fact we found no regular 
motion for velocities sufficiently below the ion acoustic velocities and for effective laser 
frequencies sufficiently below the average plasma frequency. In those chaotic situations 
the final asymptotic behavior looks like the one corresponding to uncoupled waves in the 
vacuum. The intrinsic instability of the system for too low values of the effective laser wave 
frequency, in particular, seems to preclude the formation of some special soliton solutions 
as we shall see. Chaos is present because our system is two degrees of freedom, and the 
system is two degrees of freedom because we allow for charge separation in the ion acoustic 
dynamics. Had we used quasineutrality assumptions, the system would be one degree of 



freedom and chaotic solutions would be absent. We finally point out that the transition to 
chaos we detect here is of a conservative character, so the chaotic dynamics we expect to 
see is of different type from the one present in dissipative systems P,p!2[]. 

The paper is organized as follows: in §11 we review the basic theory governing the 
interaction of a laser and an ion acoustic mode; in §111 we analyze the nonintegrable dynamics 
with help of Poincare maps and in §IV we conclude the work. 



II. THE BASIC THEORY 

A. Introductory Remarks 

We consider here the interaction of a laser wave and an ion acoustic mode in a globally 
neutral plasma consisting of light electrons and massive ions. Let us re-derive the governing 
equation for the laser mode in a slightly different way from the one used in Ref. 0. If one 
assumes that the waves are plane waves propagating along the z axis, the equation governing 
the high-frequency laser dynamics can be written in the form 

dV dz^ me M + |^|2 

where —q is the electron charge, me its rest mass, no the average density, drif. the fluctuations 
of the electronic density due to the action of the waves, and c the velocity of light. In Eq. 
(0) the laser intensity is considered strong enough to drive electrons to relativistic velocities. 
The field '^{z, t) is defined in terms of the relation 

^ = ^(x-zy)^ + c.c., (2) 

where A is the vector potential of a circularly polarized laser field, with i"^ = —1 and c.c. 
designating complex conjugate. Now we assume solutions of the form 

^{z,t) ='^{^)e'^'''-'^'^ (3) 

with k and u respectively as the wave vector and effective frequency of the laser, with 
^ = z — Vt, and ip as a real slowly varying variable; V denotes the propagation velocity. 



Next we substitute relation (j^) into Eq. (|l|) and separates the resulting equation into its 
real and imaginary parts. The imaginary component yields the relation 



from which one can determine the propagation velocity, given the frequency and wave vector 
of the carrier. The real part, on the other hand, yields a governing equation for the real 
amplitude ip, valid when V c; 



In the equation above, space is normalized by the electronic Debye length X^e = 
{Te/ATcUoq'^Y^^, time by the ion plasma frequency cu^- = ATinoq^/mi and the density by 
Ho] rrii is the ion mass, Tg is the electron temperature, /3 = c^/v^^ ^ \, /S. = uj^juj^^^ 
= Te/mei and cu^g = Airnoq'^ /rrie. Note that in our calculations, the character of the 
electromagnetic wave is determined by whether A is larger or smaller than unity. In the 
former case one has propagation in an under dense plasma and in the latter case one has 
propagation in an over dense plasma. Since it can be seen in Eq. that the coefficient 
of the linear ip-term is given by 1 — A, wave localization in over dense situations occurs 
when the natural frequency of the associated oscillator can be interpreted as an imaginary 
number. Our factor A is defined in terms of the effective frequency uj which incorporates all 
possible frequency shifts. Therefore this factor A replaces in a simplified way the slightly 
more complicated total frequency shift introduced in Ref. 0. This is why we adopt the 
present formalism. 

As for the ion acoustic field excited by ponderomotive effects associated with the electro- 
magnetic wave, we simply write down the governing equation obtained from the analysis of 
the low frequency dynamics involving the ion acoustic electric field potential the massless 
warm electronic fiuid, and the massive cold ionic fiuid. One has ||^ 



V 



c 



(4) 




(5) 



where $ has been normahzed by T^/q {q is the ion charge) and where we introduce the Mach 
number M = V/Cg with the ion acoustic velocity Cs written in the form Cg = {T^./rriiY^'^ . 
The ions are considered nonrelativistic due to their large mass. We point out here that the 
total ion and electron densities are respectively measured by the absolute values of the first 
and second term on the right hand side of Eq. (]^); in particular, the fluctuating electron 
density obtained under the assumption of massless electrons, 



5n, = e*+^-/^v^ - 1, (7) 

is the expression to be used in Eq. (||). 

Equations (|^), (^) along with (|^) govern the nonlinearly coupled dynamics of the dy- 
namical variables ip{C) and $(0- We assume propagation at a constant velocity M and take 
into consideration full nonlinear effects both in the ion acoustic and transverse relativistic 
dynamics of electrons. It is worth mentioning that alternative approaches do not restrict 
the spacetime dependence of the solutions but use weak nonlinear expansions instead p3| - 



Now, as shown in the paper by Rao et al. 0, both equations (|) and can be obtained 
from a generalized two degrees of freedom Hamiltonian H: 

H = fl-^- ^ MVM2-2$ + e*+^-^v^, (8) 

where a misprint in has been corrected. In the Hamiltonian, and P$ are the canonical 
momenta conjugate to the corresponding subscript coordinates. Hence, our system contains 
four dimensions: ipy -^V' ^ P'^- 

Since the Hamiltonian does not depend explicitly on the "time" coordinate S,-, it is a 
constant of motion. We will be interested in finding solitary waves and all those solutions 
towards which initial solitary waves can evolve in time if they are unstable. Therefore we 
shall work on the particular energy surface which allows for the presence of the configuration 
$ = '?/; = P$ = = 0, since this configuration is in fact the appropriate asymptotic solution 
for solitary pulses. We recall that from the canonical equations generated by H , Eq. (H), it 
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follows d^/d^ = — P$ and dip/d^ = P^/f3'^, and that the boundary conditions for a solitary 
pulse are, in the more traditional form, ip, d^/d^, dip/dC, — > as |.^| oo. The above 
considerations enable one to determine the constant numerical value of the Hamiltonian; it 
reads H = 1 + M^. 

The Hamiltonian supports pure ion acoustic waves; it is easy to see that regardless the 
values of P$ and $, if = 0) = P^{^ = 0) = then ipi^) = P^{^) = at any ^. This is 
not true for electromagnetic waves. Even if one starts with a laser pulse in the absence of 
any space charge fields ($ = P$ = 0), Eq. indicates that the dynamics of the electric 
field potential is subsequently driven by a tp"^ factor. 

We shall restrict the present analysis to subsonic cases, where M < 1. This causes the 
coefficient of the linear $-term in Eq.(|^) to assume only negative values. As for the laser 
field, we focus attention on over dense plasmas with A < 1. Since the linear coefficient of 
Eq. is given by a factor (1 — A)//5, as already mentioned, the system is linearly unstable 
against the generation and propagation of electromagnetic modes. Nonlinear saturation of 
the unstable process may be responsible for the creation of solitary pulses. 



B. Adiabatic Approximation 

It has been seen that laser localization occurs in over dense plasmas where A < 1. Since 
in general A is only slightly different from the unity and since the small factor 1 — A is yet 
to be divided by the large factor j3 to obtain the coefficient of the linear term of Eq. (^, 
the conclusion is that in general the following relation may hold: 



'1 - A' 



1 

M2 



(9) 



V P 

But if such a relation does hold, it is likely that the dynamics on the ip, phase plane 
tends to be much slower than the corresponding dynamics developing on the $, P$ phase 
plane. In this limiting case one could be tempted to use the results of the center manifold 
and adiabatic theorems which say that the integration on the latter phase plane could 
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be done simply by taking Eq. (||) with ifj'^ considered as a constant factor. In addition, as 
%l) slowly evolves as a function of ^, the ^-dynamics evolves in such a way as to conserve the 
action integral (l/27r) In Fig. (|I]) we take fixed values of ip to plot contour levels 

of the driven ion acoustic Hamiltonian which is here defined as: 

H^ = -^ + MVM^ - 2 $ + e*+^-^v^. (10) 

In the adiabatic regime this is the Hamiltonian controlling the fast motion on the $, P$ 
plane. For vanishingly small values of ip'^ one can see in Fig. (|l]a) an elliptic fixed point 
at the origin and a hyperbolic point located at $ < 0; trajectories move counter-clockwise 
and the typical parameters (3 = 100 M = 0.9 are considered. Now as ip"^ grows the elliptic 
point moves toward the hyperbolic point, as seen in Fig. (|l|b). For large enough values of ip 
the elliptic and hyperbolic fixed points coalesce in an inverse saddle-node bifurcation [1T7| ; 
we shall refer to this coalescence as collapse. For completeness we display in Fig. (|I|c) the 
tp, phase space obtained from the full Hamiltonian when we set $ — * and P$ —>■ 0. 
Note that when A < 1, as in Fig. (|I]c) where we consider A = 0.98, the fixed point at the 
origin is hyperbolic; for A > 1 the fixed point would be elliptic. 

A delicate point about the adiabatic approach is that the adiabatic theory may be ex- 
pected to work relatively well only if the maximum value of ip"^ throughout the entire dy- 
namics is much smaller than the critical value for which the collapse does occur. Indeed, 
if this is the case the rotational frequency around the elliptic point on the $, plane can 
be expected to be larger enough than the characteristic time scale of the ■^/'-dynamics. But 
on the other hand if the maximum value of ip"^ becomes too close to the critical value, the 
rotational frequency tends to diminish and attain values comparable to the ip-time scales. 
The rotational frequency actually vanishes at the collapse. From this point of view, Eq. 
may not be sufficient to guarantee adiabaticity, since it is derived with basis on linearization 
procedures, where the fields are assumed to be much smaller than the maximum values they 
can actually attain as time elapses. With that in mind we now proceed to derive a validity 
condition for adiabaticity that takes into account the nonlinear effects associated with the 
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finiteness of V'^- We shall see that the adiabatic range is in fact much smaller than the one 
suggested by Eq. (|^). 

Let us then examine the adiabatic trajectory of the elliptic point on the $, P$ plane. Our 
interest lies on the fact that the existence condition for this point provides a reasonably good 
estimate on the range of validity of the adiabatic regime; we emphasize that adiabaticity is 
expected to break down when the elliptic and hyperbolic points cease to exist. One can use 
the results from the center manifold theorem to estimate the location of slowly moving fixed 
points, 



\ fixed OD Ifix^d, 



9$ dP^ 

For a given value of ip, one thus have 



0. (11) 



from which a series expansion yields a relation correct up to quadratic terms 
,2/ 1\ 1-2/3-2 M2 + M^, 

pv'Jpj + W^f' '''''' ' ^ ^ 

as for P$ one has P<!>jixed = 0. We point out that in the extreme adiabatic limit where 
frequency shifts are vanishingly small, A — ^ 1, the first term on the right hand side of Eq. 
(|13|) coincides with expression (35) of Ref. 

Given ip, the quadratic relation above furnishes two roots in the variable $ if the appro- 
priate discriminant is positive. One of the roots represents the elliptic point, we shall call 
it ^eiiii^"^), and the other represents the hyperbolic point, ^hypii^^), both seen in Figs. (|l]a) 
and (03). What must be done now is to substitute the adiabatic relation ^eii = ^eiii'^P'^) into 
Eq. (^ to determine and examine the slow dynamics on the ip, phase plane. Assuming 
for a moment |$| and tp'^ small, which shall be seen to be true if A is sufficiently close to 
the unity, we drop the quadratic $-term in Eq. (|T3|) and obtain 



which shows that $ < if M < 1. Therefore in the present approximation the ^-dynamics 
is commanded by the following effective potential written, apart from a global multiplicative 

constant, as: 

VeffW = l{A-i)r + l{i+f3- r. (15) 

It is thus seen that near the ion acoustic resonance where M ^ 1, the electric field potential 
response satisfies |$| ^ t/'^ and essentially determines the adiabatically saturated value of 
the laser field: 

2 (A - 1) (M2 - 1) 

(16) 

The approximation we have used in Eqs. (|15D and (0) assumes that the quadratic $-term 



in relation (|T3D is much smaller than the others. This is true only if one is sufficiently away 
from that situation where the discriminant vanishes causing the collapsing of elliptic and 
hyperbolic points. Near the collapse, in particular, the adiabatic approximation is expected 
to break down. Given M and (3, use of Eq. (|r^) and the self consistency requirement of a 
nonnegative discriminant for relation (|T3|) finally yields a complicated relation that can be 
used as an estimate for the critical value of A where the approximation (0) is no longer 
valid 

= 4 (1-2/5-2M2 + M4) ^^^^ 

As A starts to get too close to Acr the adiabatic approximations are expected to get poorer 
and poorer. Condition ([T7|) is a rough estimate which could be refined with more detailed 
algebraic work. However we shall take it as sufficiently accurate and complete for our 
purposes. In any case the ultimate answer is yet to be given by numerical work as we will 
do next. 

III. TRANSITION FROM ADIABATIC TO CHAOTIC REGIMES 
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A. Testing Adiabaticity 



In all the following numerical applications we use /3 = 100 and M = 0.9 which yields 
~ 0.976. We promptly conclude that the validity range for the adiabatic regime is quite 
narrow, as a matter of fact, much narrower than the range predicted by Eq. (j^). Indeed, 
for the chosen values of (3 and M, and considering over dense plasmas, Eq. (^) basically 
imposes no essential restriction on the value of A, a failure of the linear theory as mentioned 
before. In the following we shall see that the estimate based on Acr, Eq. (|r^), is much more 
accurate than the one based on Eq. (^, and that the destruction of the adiabatic regime is 
in fact associated with a transition to chaos. 

Before examining the validity ranges and the transition, let us first perform some initial 
simulations of Eqs. @ and @ to make sure that the adiabatic regime is in fact present if 
condition (|l^) is safely observed. To do so we start a single initial condition with A = 0.99, 
and with 7/; = P$ = 0, $ = 0.0005 and = 0.0024 such that H = l + M"^. We plot the time 
series for ip{C,) and $(^) in Fig. (Q). In the figure we see that $ undergoes a fast oscillatory 
motion while ip evolves in a much slower time scale. The adiabatic features can be also 
visualized on the ip,Pip phase plane as in Fig. (^. In the figure we compare three solitary 
trajectories: the exact trajectory, the adiabatic trajectory calculated from Eq. under the 
assumption (0), and the trajectory calculated from Eq. (j^) under the assumption $ 0. 
The adiabatic trajectory yields a fairly good approximation to the actual trajectory. Here 
we use A = 0.98. This chosen value of A is slightly smaller than in the previous figure, 
because it allows a clearer view of the differences between adiabatic and exact trajectories. 
For values of A closer to the unity the approximation gets better and better until such a 
point where no distinction can be appreciated. 
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B. The Transition 



Now the question refers to what happens as the parameters are varied beyond the vahdity 
range for the adiabatic regime. To simphfy the discussion we shall focus attention on the 
behavior of the system as A decreases. As our system is Hamiltonian (see Eq. (^) with 
two degrees of freedom, we make use of the Poincare map methodology and plot the pair of 
phase variables $ and Pq> each time = with dP^/dC, > 0. Several initial conditions are 
launched with the numerical values for = 0) distributed within a small range typically 
satisfying —0.001 < = 0) < 0.001. Similarly to the initial simulation presented before, 
for all initial conditions we always take ip = P<i, = and calculate the corresponding initial 
P^ from the constant numerical value of the Hamiltonian, H = 1 + M^. This kind of 
launching conditions initially places the system in the vicinity of the solitary solution which 
is the solution containing the point ^ = P^ = ip = P^p = 0. In integrable cases the ensuing 
nearly solitary trajectories progress in fact as trains of solitons, but even in the nonintegrable 
cases where solitons are not seen, the trajectories still cross the P^ = plane several times, 
an essential condition for the construction of the maps. 

We start by displaying in Fig. (|a) the map obtained when A = 0.98. For such a value 
of A the adiabatic regime is expected to prevail. In agreement with that, what is seen in 
the plot is a set of regular orbital concentric curves. Note that the elliptic point appears 
to be located at a negative value of $ simply because this is the value of the electric field 
potential when the recording conditions P^ = 0, dP^/d^ > are satisfied. Now if one starts 
to decrease A the transition to chaos is expected to occur. Let us move on to Fig. (^b) where 
A = 0.975. As anticipated from the analytical estimates, a considerable amount of chaotic 
activity can already be identified. The central region of the map is completely surrounded 
by a blend of stochastic orbits and resonant islands. In particular, it appears that the soliton 
solution which corresponds to the central fixed point no longer exists. In Fig. (^) we enlarge 
part of Fig. (^3) to show details of the resonant islands. These small remaining regions of 
regularity of the phase space are then totally suppressed when one reduces A further below. 
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In Fig. (§d), for instance, we consider the case A = 0.97 to show a deep chaotic regime. 
In conclusion, a complete destruction of the regular trapping region does indeed occur as A 
decreases. But if there is no trapping region, how would the trajectories behave? We will 
address this issue next. 



C. Asymptotic States Beyond the Transition 

We can actually see some persistency in chaos only for those intermediary situations 
where A is not too close to the unity, but not too small. Indeed, if A becomes sufficiently 
small, say A < Acr, no trapping region is effectively formed on the $, P$ plane. Even initial 
conditions originally launched within the trapping region predicted by the adiabatic theory 
when ip = 0, do not remain there. As a matter of fact, the trajectories are eventually ejected 
into unbounded regions of the phase space when A crosses the validity limits of adiabaticity. 
For those cases, chaos would be at most a transient which would take place during initial 
instants, before ejection. The question to be asked now should be on the configuration of 
these unbounded orbits. What we have observed is that once the trajectory escapes from 
the trapping region in the ion acoustic phase space, it starts to follow the open fiow lines 
of Fig. (P. This is confirmed in Fig. (^) where we show a continuous plot displaying a 
freed trajectory which was started with A = 0.97 and = 0) = 0.0001 (figures (^) and 
(Pd) are not Poincare maps; the trajectory points are periodically recorded with a small but 
constant time step). The fact that the subsequent trajectory evolves along the fiow lines 
implies that $ gets more and more negative. Now, if one considers Eq. (|^) one readily sees 
that regardless the value of ip, 

lim H = ^ + _El + M VM2 - 2 $. (18) 

In other words, laser and ion acoustic fields become decoupled in this limit. Since AjS > 
then the corresponding dynamics of the laser field must necessarily become that of an 
undriven harmonic oscillator. This is what is shown in Fig. (^) where we project the 
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same dynamics of Fig. (^) now on the ip, phase plane. After a certain amount of 
time following a figure eight shape like the ones seen in Figs. ([I|c) and there is a 
dynamical transition to the circular shape so characteristic of the harmonic oscillator. The 
instant of the transition coincides, as it should, with the moment of ejection seen in Fig. 
(^a), and occurs approximately after 30 cycles of the laser wave in its initial figure eight 
phase trajectory. Another interesting point connected to this asymptotic state is that as 
$ ^ — oo, the particle density becomes very small (see Eq. (0)). Noticing that $ < 0, the 
ponderomotive field created by the laser induces an initial potential well in which interior 
the ion fiuid undergoes acceleration becoming less dense. If the amplitude of the laser is too 
large the process is unstable and never arrests. When the density is low the laser becomes 
an almost standing wave with a small propagation velocity of the crests {V <^ c). 



D. The Role of Relativistic Effects 

Our original equation, Eq. (^), includes full relativistic electronic nonlinearities. As a 
final topic it is perhaps interesting to discuss the role of these relativistic nonlinear effects 
as compared to ponderomotive nonlinear effects. 

What we find here is that while in the adiabatic regime, saturation is essentially gov- 
erned by ponderomotive nonlinearities. This is the basic conclusion associated with Eq. 
(|T6|). Now we would like to know whether or not relativistic effects grow in importance 
in chaotic regimes. To this end we perform two pair of simulations that are displayed in 
Fig. d^). In the upper panel of each pair we depict time series for considering exact 
fully relativistic nonlinear dynamics. In the lower panel we consider the time series with 
relativistic mass correction suppressed. It is seen that in the adiabatic regime of Fig. (||a) 
where we consider A = 0.98, relativistic effects are not prominent, as both figures are almost 
identical. On the other hand, in Fig. (^o) where we consider A = 0.97 relativistic effects 
are of relevance. Indeed, in the time series of Fig. (^) we see that for the particular initial 
conditions we use, the nonrelativistic trajectory undergoes much earlier ejection from the 
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initial trapping region. For other initial conditions, the escape order may be reversed and the 
nonrelativistic trajectory may be ejected after the relativistic one. What is really remarkable 
here, however, is that even in our case of small field amplitudes where relativistic effects are 
small, relativistic and nonrelativistic trajectories may largely differ in chaotic regimes, the 
reason for this being the extreme sensibility of chaotic systems to parameter variations. 

IV. FINAL CONCLUSIONS 

We have performed a nonlinear analysis on the interaction of high power laser waves with 
ion acoustic waves in a plasma. We assume stationary propagation in an over dense plasma 
and consequently show that three generic configurations take place. If the effective laser 
frequency is only slightly lower than the average electron plasma frequency, the ensuing 
dynamics is adiabatic. The faster varying ion acoustic field is adiabatically enslaved to 
the slowly varying envelope laser field, and the resulting electromagnetic envelope solitons 
are likely to exist and remain stable. As one starts to decrease the effective laser frequency, 
adiabaticity becomes progressively poorer. A blend of confined chaotic regions and nonlinear 
resonance islands are seen on the $, P$ phase space. Then, for yet smaller values of the laser 
frequency adiabaticity is completely destroyed. Initial conditions are rapidly ejected from 
the trapping region on the ion acoustic phase space and proceed to move along unbounded 
curves towards $ —>■ — oo. In this limit the system becomes decoupled, the laser field starts 
to behave like a vacuum field and the plasma becomes progressively rarefied. 

Relativistic effects are moderate in adiabatic regimes but of considerable relevance in 
chaotic regimes. We saw that trajectories with relativistic corrections artificially removed 
can largely differ from the exact ones, the reason being the sensibility of chaotic systems to 
parameter variations. 

One general conclusion obtained here is that the parameter range for the existence of 
solitary wave is extremely narrow. This numerically confirmed fact can be predicted only 
by nonlinear estimates like that provided by Eq. (]T7|), which takes into account the effect 
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of finite values of if)'^ on the existence of fixed points. Linear estimates like that of Eq. 
(y) are much less accurate. Due to this narrow existence range, we have not observed, for 
instance, double hump solitons like those obtained in Ref. 0. The point is that according 
to the calculations done in this reference, and making the appropriate connections and 
translations between the various formulas, double hump solitons exist only when A reaches 
small values, A ~ 0.7. However, for such a small values of this parameter, our system has 
already lost stability due to the transition to chaos. 
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FIGURES 

FIG. 1. Phase space on the P$ plane for fixed values of if), with M = 0.9 and /3 = 100; ijj = 

in (a), and ip = 0.001 in (b). In (c) we display the tp,P.,p phase space using A = 0.98. 



FIG. 2. Time series in the adiabatic regime; M = 0.9, P = 100, A = 0.99 

FIG. 3. Accuracy of the adiabatic approximation tested on the ip, plane; M = 0.9, /3 = 100 
and A = 0.98. 

FIG. 4. Transition to chaos for M = 0.9 and P = 100; A = 0.98 in (a), A = 0.975 in (b), and 
A = 0.97 in (d). In (c) we show details of the resonant islands seen in (b). 

FIG. 5. Asymptotic states on the phase planes. The escaping trajectory seen in (a) and the 
circular trajectory seen in (b) do appear on the respective phase planes after approximately 10 
figure eight cycles of the i/j, variables; M = 0.9, /3 = 100, A = 0.97 

FIG. 6. Influence of relativistic nonlinear ities on the dynamics. In the lower panel of each 
figure relativistic effect is artificially suppressed. In (a) we consider an adiabatic regular regime 
with A = 0.98 and in (b) we take A = 0.97 to produce an escaping trajectory after some transient 
chaos. It is seen that in the chaotic case one should not discard relativistic effects. 
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Fig. 4b: Rizzato, Lopes, and Chian 
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